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Abstract 

We study coupled maps on a Cayley tree, with local (nearest-neighbor) 
interactions, and with a variety of boundary conditions. The homogeneous 
state (where every lattice site has the same value) and the node-synchronized 
state (where sites of a given generation have the same value) are both shown to 
occur for particular values of the parameters and coupling constants. We study 
the stability of these states and their domains of attraction. As the number of 
sites that become synchronized is much higher compared to that on a regular 
lattice, control is easier to effect. A general procedure is given to deduce the 
eigenvalue spectrum for these states. Perturbations of the synchronized state 
lead to different spatio-temporal structures. We find that a mean-field like 
treatment is valid on this (effectively infinite dimensional) lattice. 

PACS number(s):05.45. +b, 47.20. Ky 
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1 Introduction: 



Coupled map lattices (CML) have been explored in a variety of contexts in recent 
years, particularly as prototypes of spatially extended systems. These are simple 
models wherein both space and time play a role; furthermore, it is anticipated that 
the insight gained over the past two decades in studying low-dimensional nonlinear 
dynamical systems can be profitably exploited in providing an understanding of such 
complex systems]]]]. 

The phenomenology displayed by coupled maps on regular one and two dimen- 
sional lattices has been extensively studied by Kaneko|J. In addition, CML have been 
used to model a wide variety of complex phenomena, such as the study of the kinetics 
of phase ordering processes 0, crystal growth ||, neuronal systems ||, optical fibres 
H and pattern formation ]7]]. Chate and Manneville have also used CML to model 
spatiotemporal intermittency||. A route to a spatiotemporally inhomogeneous state 
through wavelength doubling bifurcations has also been recently identified]^]. 

In this paper, we study coupled maps on a Cayley tree. This lattice is embedded 
in infinite dimensions, and thus should give some indication of CML phenomenology 



in higher dimensions ]10[. Although the Cayley tree (sometimes termed the Bethe 
lattice) is an idealized hierarchical lattice with no immediate physical application, 
it is convenient for study since there are no closed loops. Furthermore, the Bethe 
lattice is the simplest sort of branching media model encountered in many physical 
processes. 

Previous studies of coupled map systems have largely (except for a study by 
Cosenza and Kapral]ll| of CML on a Sierpinski gasket) been carried on with local 
coupling on regular lattices in 1- and 2- dimensions, or with global coupling, in which 
case there is no notion of lattice geometry. Our motivation in choosing the Bethe lat- 
tice is two-fold. Apart from the mathematical convenience, it is worth considering that 
in many physical situations, the medium supporting dynamics could be nonuniform; 



in cases like chemical reactions in porous media or on diffusion-limited-aggregation 
clusters, heterogeneity can lead to a hierarchical structures]!^]. We note that hier- 
archical structures have long since been studied also in spatiotemporal systems like 
neural nets also due to their exponentially higher storage capacity JT3]]. 

A related question of some current interest is the control of macroscopically cas- 
caded dynamical systems. The synchronization of a large set of oscillators connected 
in series [ 



14, [13, [L6L O] i n a given geometry and with particular boundary conditions, 



is directly related to the problem of whether a similar CML can support a synchro- 
nized state[T8|. We address this problem, and show below that the criterion in CML 
for a synchronised (but chaotic) state to be stable is that only one Lyapunov exponent 
is positive and all the rest are negative. In both the cases, i.e. synchronization of the 
coupled oscillators or coupled maps the essence of the problem lies in the nature of the 
eigenvalues and eigenvectors of the interaction matrix. In the present work we deal 
with the situation of asymmetric coupling which is easily obtained in experiments [T"5| . 

The plan of this article is as follows. We define our model, and the boundary 
conditions, and show that the stability of synchronized states depends on the spectrum 
of eigenvalues of the interaction matrix. This is related to the connectivity matrix 
for the lattice, and has a singular-continuous structure. Different patterns can arise 
from the secondary instabilities. We note that for coupled piecewise linear maps the 
study of the interaction matrix gives the whole Lyapunov spectrum, which is related 
to the correlation length. 

In several recent studies of globally coupled maps[0, a breakdown of the "law 
of large numbers" has been observed. We find that for this system, even with local 
coupling, the mean-field description is valid in the macroscopic limit. This is probably 
related to the infinite dimensional nature of the lattice. 
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2 The Model: 



We study CML on a Cayley tree with coordination number 3, which we construct in 
stages as follows. At the first level, there are three branches, each of which split in 
two and each of them further bifurcates and so on. Each site thus has one parent 
and two daughters, except for the origin which has no parent site and 3 daughters, 
and boundary sites which have no daughters. Each site on the lattice is assigned 
a variable x, which evolves according to a deterministic rule depending on its own 
value and the value of the nearest neighbors. The evolution rule is taken to be the 
following: 

x(i, t + l) = h f(x(i, t)) + h p f(x(i p , t)) + h d Y, fW, t)) (1) 

where / is the function that determines how the lattice variables evolve. A common 
choice for / is the logistic map, f(x) = /ix(l — x). The notation above is as follows: i p 
is the parent of site i, while id are the daughters. h , h p and h d are coupling constants 
and are taken such that h + h p + 2hd = 1. Thus the evolution is contained in the 
same phase space as that of a single map, e.g. [0,1] in case of logistic map. We also 
assume that the couplings are positive though most of our results do not explicitly 
require this. 

The evolution rules for the origin and for boundary sites are somewhat different, 
since the former does not have a parent and the latter do not have daughters. For 
the origin, 

x(0, t+l) = o f(x(0, t)) + o' £ f«J, t)) (2) 

jeo d 

while for the boundary sites, 

x(k, t+l) = b f(x(k, t)) + b'f(x(k p , t)) (3) 

with different choices for oq, bo, o', b' to ensure that the dynamics remains in the same 
phase space. 
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Each point on the Bethe lattice can be indexed by a string a = {ai, 02, ... , 
where % denotes generation. a\ can take the single value a\ = 0, a<i can be 0, 1 or 2, 
and for j > 2, the a/s take two values, a 3 ■, — or 1. With this notation, it is possible 
to assign a unique number n(A) = n(ai, . . . , a*) to each lattice point. 

n(ai, a 2 , . . • , a*) = #(i - 1) + a* + 2a;_i + 2 2 a;_ 2 + . . . + 2 i ~ 2 a 2 + 1 (4) 

where 

^(<) = 1 + 3C2*- 1 - 1) (5) 

is the total number of sites at the ith generation and g(0)=0. (See Fig. 1). We can 
thus formally write 

X(t + 1) = IF(X(t)) (6) 

where F(X(t)) is a column matrix, X(t) is an array of variable values assigned to the 
lattice points arranged in ascending order of n(A), and I is the interaction matrix. 
For example, for the Cayley tree with 3 generations i. e. 10 sites, the equation above 
reduces to 



X(t + 1) = 
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(7) 



where 



(8) 



F(X(t)) = (f(x(l,t)), f(x(2, t)), f(x(3, t)),..., f(x(g(k),t)Y 

where T denotes transpose. Note that the interaction matrix / is such that the 
evolution can be written as 



S(M+I)=£jy/(*(J\0) 



(9) 



where the sites have been labeled by the unique number n(A). 



5 



3 Coherent Patterns and Their Stability: 

Inspection of the symmetries of the lattice allow for the determination of allowed 
coherent patterns. It is easy to verify that if one starts with the pattern in which 
all the points at each generation % are assigned the same value Zi(t) at time t, the 
nature of the pattern cannot change in time since points at the same generation have 
equivalent evolution rules and remain synchronized. With parameters 

o + 3o' = b + b' = h + h p + 2h d = 1 (10) 

another simple pattern is possible. This is also node-homogeneous, with Zi = z for all 
i, i.e. all the points on the lattice are synchronized since the evolution is essentially 
that of a single map /. 

These 'allowed' patterns can be observed in practice if and only if they are linearly 
and convectively stable against small perturbations. In the present work, we have 
mainly dealt with linear stability analysis of this system in the stationary frame, and 
while we have not analysed convective stability, numerical experiments suggest that 
no extra instabilities other than the ones in an equivalent one dimensional model 
creep in. This directly evolves on the fact that there are no loops on the lattice- 
there is only one direction in which instabilities can be enhanced in a moving frame 
of reference, and these are the same as in the equivalent 1-d model. 

For the linear stability analysis the eigenvalues and eigenvectors of the matrix J = 
lim^oo J(r) where J(r) = J T . . . J 2 J±, are (asymptotically) relevant. The Jacobian 
matrix at time t, i.e. J t is given by Jt(i,j) = I(i,j)f'(xj(t)) and Xj(t) = x(t) for all j. 
Thus the Jacobian matrix is J = lim T ^ OQ [I] T f \x T ) f '(x T _i) . . . f'(xi). The eigenval- 
ues of J are lim^oo A[ where Aj = Vi\ where Vi, % = 1,2, ... , g{k) are the eigenvalues 
of the interaction matrix I and A = lim^oo \f'(x(t))f'(x(t — 1)) . . . f'{x{l))\ l / t . The 
relevant eigenvectors are those of /, and the problem reduces to a study of the eigen- 
values and eigenvectors of the interaction matrix. 
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The fact that coherent patterns are allowed (by the condition in Eq. |1(]) implies 
that a right eigenvector of the interaction matrix is e\ = [1, 1, . . . , 1]. This is a 
characteristic of row stochastic matrices, and corresponds to the eigenvalue A for the 
product of the J's. From Greshgorin's theorem|22| this is the largest eigenvalue. 
Consider a small deviation, Ao = [8±, 8 2 , ■ ■ ■ , S g (k)] from the homogeneous pattern 
[x, x, . . .]. We can reexpress A in the basis of eigenvectors e±, e 2 , . . . , e g ^) 

A = a x ei + a 2 e 2 + ■■■ + a g ( k )e g ( k ). (11) 

After t iterations the deviation from the homogeneous condition will be 

A t = diA^ei + a 2 \ t 2 e2 + . . . + a g {k)X g (k) e g{k)- ( 12 ) 

If the only eigenvalue with modulus greater than unity is Ai = A and \\Vj\ < 1 for 
j > 1, i.e. the rest are less than unity in magnitude, then for large enough t we can 
write 

A* ~ aiA* ei (13) 

The perturbation grows along the direction t\ = [1,1,..., 1] and any random devia- 
tion will eventually be homogenized. 

Thus the necessary (though not sufficient) condition for the synchronized pattern 
to exist (and evolve chaotically in time) is that Ai is the only eigenvalue greater 
than unity and all others are less than unity in magnitude: a linearly stable coherent 
pattern- in the infinite lattice limit- therefore requires a finite gap in the eigenvalue 
spectrum of the interaction matrix. 

The interaction matrix is analogous to the tight-binding Hamiltonian on the Bethe 



lattice pi I although the eigenvectors are different (since the matrix is not necessarily 



symmetric or hermitian). However, using similar arguments||21||, it can be seen that 
all the sites at a given generation are equivalent in the sense that if sites at every 
generation are synchronised, this pattern will continue to exist in absence of small 
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perturbations or noise, since the evolution rule is the same for all of them. Further- 
more, one can see that if any two sites which have the same parent are interchanged 
along with their sub-trees, the system is left unchanged. Using the equivalence of all 
points at a generation and the permutation symmetries of the lattice, the similarity 
transformation that will block-diagonalize the interaction matrix can be deduced to 
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(The first three vectors follow from the fact that lattice points at each generation 
are equivalent. The fourth and sixth vectors simply represent the two linearly inde- 
pendent and mutually orthogonal interchanges possible between points at the second 
generation, ((1, 0, —1) ± (0, 1, —1)) while the fifth and seventh vectors are similar in- 
terchanges within siblings with the phase derived from the parent site. The last three 
vectors arise from the interchange between the siblings of the same parent.) 

Thus the block diagonalizing matrix is written using permutation symmetries of 
the underlying lattice; the blocks are as follows. 



The two doubly degenerate eigenvalues are the eigenvalues of the matrix given below. 
They correspond to the fact that one can have two independent permutations in the 
three branches at the first node. 



be 




(15) 




(16) 
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Finally we have a triply degenerate eigenvalue is bo (reflecting the fact that one can 
have permutations between the daughters of any of the three branches at the second 
node without affecting the matrix). One can see that the consecutive blocks giving 
eigenvalues are just like earlier blocks except that first row and column of the matrix 
are removed. This construction can be trivially extended to a matrix of higher order. 
The matrix S^LS is block diagonal. 

This scheme can be generalized to higher dimensions and the diagonalizing matrix 
for the kth stage can be deduced as following. Specifying the non-zero components 
of the column vectors (in the notation of Eq. [| to denote the components) the first k 
vectors are as follows, 

V n(ai) = 1> 

V l(a u a 2 ) = 1, (17) 



. . . and 



J n(a 1 ,a 2 ) 
u n(ai,a 2 ,a 3 ) x ' 



1) k = 1 (1 i 

<y n(ai,a2,a 3 ...,a fc ) ■ l xc 



{e.g. the first 3 columns of the matrix defined in Eq. H) Then we have two sets of 
k — 1 vectors. One is 



U,k+1 _ k+1 _ -i k+l _ _ 2 l 
\ V n(0,l) ~ U n(0,2) — L > ^(0,3) ~~ Z J ' 
( v k+2 _ k+2 _ 1 k+2 _ _ 2 | 



(19) 



. . . and 



r k+{k-l) _ k+(k-l) _ 1 k+(k-l) z™-* 

\ V n(0,l,a 3 ,...a k ) ~ V n(0,2,a 3 ,...a k ) ~ X ' V n(0,3,a 3 ,...a k ) ~ Z f- l ZU J 



(See e.g. the fourth and fifth columns in Eq. [14].) The other set is 

{^(0,1) = ^i V n k (0,2) = ~ 1}> 
| ? .2fc+l _ 1 2fc+l _ _1\ 
l U n(0,l,a 3 ) ~~ U n(0,2,a 3 ) ~ l 5 



(21) 



. . . and 

%i(0,l,a 3 ,...,a fc ) — A -^ v n{0,2,a 3 ,...,a k ) 



r 2*-l+(fc-l) _ 2 fe-l+(fe-l) _ n (22) 



(See e.g. the sixth and seventh columns in Eq. [TJ].) The next three sets of k — 2 



vectors each are given as follows. The first is of the type 

{^(0,1,1) = l> V n(0,l,2) = — 1}' /.-,■-)., 
f v 3k _ 1 3fc — _U ^ > 

1^(0,1,1,04) ~~ L i ^(0,1,2,04) ~~ J ' 
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and 



{v 



The second is 



. . . and 



The last set is 
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(24) 



(25) 



(26) 



{ V n(0,3,l) ~~ 1' ^(0,3,2) — 1} 

and so on. (The last three columns in Eq. [14] are i; 3 * -1 , t> 4fe-3 and v 5k ~ 5 . For > 3 
newer sets will appear.) Now we will have sets of vectors that will give blocks of size 
A; — 3. The next six blocks of k — 3 vectors arise from permutations between the points 
on the fourth generation and their descendants and are of the same type as the three 
sets of k — 2 vectors mentioned above which result from three independent permu- 
tations possible between the six points on the third generation. One can continue 
this scheme untill reaching the boundary. The number of points on the boundary is 
g(k) — g(k — 1) = 3 2 k ~ 2 and (g{k) — g(k — l))/2 = 3 2 fe_3 permutation vectors are 
possible (see Eq. H) which will give a block of size 1 with the same degeneracy as 
the number of permutations possible on the boundary. 

For boundary conditions in which o Q = h + h p , d = 2/i rf /3, b = h + 2hd, b' = h p , 
at stage k, the first block of the block-diagonal form is 



/ ho + h p 2h 
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h n 2h n 
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ho 2h d 
... h p h + 2h d J 



(27) 



which exploits the equivalence symmetry of all the sites at a given generation. The 
second block, which exploits the permutation symmetry of the points on the first 
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generation is 

( h 2h d . . . 0\ 

h p h 2h d 

h 2h d 

\ ... h p h + 2h d J 

and so on. The last two blocks are 



(2* 



h 2h d 
hp h + 2h d 



( h + 2h d ) (29) 



The first block of order k appears once in the block diagonal form, the second of order 
k — 1 appears twice, and next blocks of order k — n(k — l>n> 2)appears 3 2™~ 2 
times. The first k eigenvalues are therefore nondegenerate, then k — 1 eigenvalues are 
doubly degenerate, k — n eigenvalues have degeneracy 3 2 n ~ 2 for k — 1 > n > 2 etc. 

From the structure of the matrices and their degeneracies, one sees that for the 
Cayley tree with one more generation, the k — 1 degenerate blocks are retained (with, 
however, the doubly degenerate block becoming triply degenerate and other blocks 
doubling their degeneracy) with an additional block which has degeneracy 2. The 
block corresponding to non-degenerate eigenvalues is, however, completely changed. 
The density of states has to be singular continuous since the new eigenvalues that 
are created have a lower degeneracy: the eigenvalue spectrum is a sum of delta peaks 



and is nowhere differentiable, as is common in hierarchical systems jy] . 

In the situation when the synchronized state is linearly stable in the stationary 
frame, the typical degeneracy structure of eigenvalues is as shown in Fig. 2, for the 
parameters and boundary conditions as discussed below. The structure is generic if 
the system is linearly stable, but the width of the gap varies with the parameters. 
For piecewise linear maps, e.g. f(x) = rx mod y, the Jacobian is constant in time, 
and the spectrum of eigenvalues of the interaction matrix determines the Lyapunov 
spectrum of the CML, and thus (via the Lyapunov dimension) the fractal dimension. 

We can see from the degeneracy structure that about a quarter of the eigenstates 
have their support fully from the boundary. The next layer is approximately half 
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of this number, and so on, with the number of states which have their support up 
to a length I from the boundary reducing exponentially. This is in keeping with the 
expectation that the rate must be faster than that in finite dimensional spaces where 
the number of modes with wavenumber \k\ < k is proportional to n d . 

Note that in the block-diagonal matrix, the blocks are tridiagonal and (for positive 
couplings) all elements are positive. Such a matrix can be transformed to symmetric 
form|22]| and thus all its eigenvalues are real: there can be no Hopf bifurcation leading 
to the instability of a synchronized state. 

Since the consecutive blocks are the principal (tridiagonal) submatrices of the 
earlier block, the eigenvalues are interlaced p3|. In other words, the bounds for the 
eigenvalues of the lower block are contained in the bounds for the eigenvalues for the 
higher block, and it is enough to consider the first two blocks in order to study the 
stability of spatially synchronized state. 

For the first block, the non-degenerate eigenvalues are given by h + h p + 2hd which 



is set to 1 by definition and the other k — 1 eigenvalues are h$ + 2^j2hdh p cos{9) where 
9 = 2ni/k, i = 1, 2, . . . , k — 1. The eigenvalues will have a gap if 2hd ^ h p . 

Consider the second block of order m — k — 1, which is tridiagonal, and can 
be symmetrized using a similarity transformation involving a diagonal matrix with 



elements D^i = [^J2hd/h p ] 1 1 . This yields a tridiagonal matrix O such that the 
diagonal elements remain unchanged and all the elements on upper and lower diagonal 



are J2hdh p . Using Greshgorin's theorem^] again, one can see that the largest 



eigenvalue cannot exceed h + 2^2h p h d if h p > 2h d . As explained above, the analysis 
of the first two blocks suffices to explore the stability of synchronized state and thus 
the other blocks do not modify the gap in the eigenvalue spectrum of the first block 



if h p > 2hd- 



Aranson, Golomb and Sompolinsky||14jj consider asymmetrically coupled 1-d chains 
with open boundary conditions where there is a convective instability of synchronized 
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patterns: perturbations grow in the moving frame of reference, destroying macroscopic 
coherence. As we have shown above, under these conditions macroscopic chaos is lin- 
early stable in stationary frame also on the Cayley tree. However, the difficulty in 
synchronizing large systems is less pronounced in this case. Due to the ultrametric 
topology, much larger systems can be synchronized under the conditions above. With 
open boundary conditions and asymmetrical coupling, coherence is more easily estab- 
lished in the present case. For example, for hi — 0.7 in one direction and h% — 0.1 in 



the other direction AransonJTj] have a coherence length of around 55 for the choice 
of map f(x) = a — x 2 with a value of a such that the eigenvalue for a single map 
is 1.26. Taking h p — .7, 2hd = .1, Fig. 2 shows a plot of eigenvalues as a function 
of degeneracies at these parameter values for 50 generations. One can clearly see 
the gap between a single nondegenerate eigenvalue above unity and the others below 
unity. For the above parameters we can easily obtain a coherent pattern for k = 20 
with random initial conditions: a CML with « 10 6 sites is easily synchronized [2~3] to 
within 10~ 4 even under single precision (16 binary digits) evolution! This is in stark 
contrast to 1-d coupled CML, which has a coherence length of about 55 sites. This 
example is a good illustration of the dramatic increase in stability with hierarchical 
connectivity. 

To check that no other instabilities than the ones expected from an equivalent one- 
dimensional model come into picture, we looked at the function f(x) = 1.39 x mod 1 
with the same choice of coupling constants as above. Here the coherence is within 
10~ 5 for first 6 sites on a lattice of 1-dimensional lattice and even on Cayley tree it is 
maintained for 6 generations. This is expected, since there are no closed loops and the 
only direction in which the instabilities can flow and grow is the one from the centre 
to the boundary. However, we can see that since the number of sites synchronized is 
exponentially higher on a tree with equivalent generations than an one- dimensional 
lattice, exponentially larger number of sites are synchronized on trees at equivalent 
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parameters. The base of exponent is related to the number of branches. 



Auerbach||15|| has shown that one can circumvent the difficulty arising from con- 
vective instabilities on a 1-d lattice by using system-size dependent feedback control. 
In essence, we achieve the same ends through a change in geometry, without extra 
controls. The boundary conditions oq = bo = ho, o p = 3bd = h p + 2h s also gives the 
same result, which indicates that some more variants are possible for open boundary 
conditions and asymmetric coupling. 

Now consider the node-homogeneous pattern. The stability matrix is given by 



j = UJ t 



(30) 



t=i 



Thus analysis of the eigenvalues of the product of matrices is reduced to the analysis 
of the eigenvalues of the product of blocks. This is a great simplification since instead 
of considering matrices of order 2 k , where k is the number of generations, we need 
consider k matrices of order k and below. The analysis of the Jacobian matrix reduces 
to analysis of the matrices 
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(31) 



•• V b ) V f'[z k (t)] J 
and so on. 

Again the degeneracy structure is the same as for the interaction matrix, and 
the Lyapunov spectrum is the sum of delta peaks and is a everywhere discontinuous 
function, as for the fully synchronized state (which is a special case of the node 
homogeneous structure). We can similarly argue that the condition for stability of 
the node homogeneous state (evolving chaotically in time) is that the first block 
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corresponding to the nondegenerate eigenvectors is the only one with eigenvalues 
of modulus greater than unity, all other blocks having eigenvalues with modulus 
less than unity. (This is since the first block corresponds to eigenvectors that have 
contribution from all the generations, and the contribution from all the points of the 
same generation is the same). 

A simple example of such stable patterns can be constructed for f(x) = rx mod 
1, with boundary conditions Oo = 0,&o = 0, d = hd,b' = h p , and parameters r = 
Y / ^3)/2,/i = 0, hd = h p = 1/3. For the Cayley tree with 5 generations i.e. 46 sites, 
it can be shown that the eigenvalues are higher than unity for the first block alone. 
Numerically, one can easily get node homogeneous patterns, starting from random 
initial conditions. Thus the possible coherent patterns are characteristically different 
from those on regular lattices and the stability analysis is also distinct [p6| . 

4 The Infinite Dimensional Character: 

Now we study the properties of this model which should reflect the fact that it is 
embedded in infinite dimensions, where a mean-field-like treatment can be expected 
to be valid. A collective variable^, ^Of h{t) is defined as 

W) = "Tlx E (32) 
9\ K ) »=i 

where g(k) is the total number of sites on the Cayley tree with k generations as noted 
above, and f(x) = fix(l — x) with h = 1 — e, h p = h d = e/3, b = o = h , b f = e, 
d = e/3 while the parameter values are e = 0.1 and = 4. The return map of this 
variable i.e. h(t + 1) vs. h(t) is a filled ellipse, whose size decreases rapidly with the 
number of generations. We conjecture that in the macroscopic limit it tends to a 
fixed point: i.e. though the evolution is chaotic for the system the collective variable 
is invariant in time. The mean square deviation of h(t) decays like 1/N where iV is 
the number of sites (see Fig. 3), quite unlike the case of globally coupled mapsffZTj], 
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where some reorganization occurs in a way such that total number of independent 
degrees of freedom is not linearly proportional to the number of sites. This is not 
totally unexpected ||2 (J|| since the values being summed are not independent random 



variables. This also means that mean field is not valid in these systems. However, 
this expectation is fulfilled for the Bethe lattice although the variables that are being 
summed are not only not independent but are also not identically distributed: the 
boundary evolves differently from the bulk and boundary effects are not negligible in 
any limit since half the points reside at the boundary. Fig. 4 shows the probability 
distribution of the central sites and the boundary for the above case and they are 
clearly different. However, the sum behaves in a way that is expected from the sum 
of iid variables. (We have verified that this behaviour holds at other values of the 
parameters) . 

The recovery of mean field in infinite dimensions is interesting. Chate and Man- 
neville [f2I| have found that the mean field like approach works better in higher di- 
mensions in spatially extended systems and connectivity plays a relatively marginal 
role. Further studies will be necessary to determine the upper critical dimension for 
this problem. We are exploring this question. 

5 Conclusions: 

In summary, in this paper we have considered a coupled map lattice in an ultrametric 
space. We show that synchronized but temporally chaotic systems can be stabilized 
more easily in this space. We present a simple method of obtaining the eigenvalue 
spectrum for the node homogeneous and spatially synchronized structures using sym- 
metries. We emphasize different properties that owe their existence to the hierarchical 
connections and to the infinite embedding dimension. 



In an ecological model, Hogg et. al. |25[ have suggested that a hierarchically 
ordered random system should be more stable compared to unstructured systems, and 
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argue that the stability should scale like log (system size) instead of (system size): 
thus exponentially larger systems should be stabilized in hierarchical organization. 
Our results are in conformity with this expectation, although the system studied in 
23] has a random branching. They have also found out that assymetric interactions 



give a higher stability. This is also expected from our analysis. It is difficult to find 
such a naturally occurring system with a clearly demarcated tree-structure. However, 
the properties we are trying to emphasize are qualitatively unchanged so long as the 
connectivity remains free of loops (and the resulting feedback). 

Our results have immediate relevance to the problem of synchronization in chaotic 
systems, which is currently evoking considerable interest fl30]. In 1-d, this problem 
has been extensively discussed [p"6] , [T7| . Our results demonstrate a novel method of 
stabilizing large systems and are of practical utility in this context. 
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Figure Captions 



Fig. 1 The Cayley tree with three generations and the labeling scheme. Arrows are 
directed towards the daughters. 

Fig. 2 Eigenvalues and their degeneracies for the synchronized state for k = 50 gen- 
erations. The parameters are h p = .7, h = .2, h s = .05 and A = 1.26. One 
can clearly see the gap that separates a nondegenerate eigenvalue greater than 
unity and all the others below unity. Degeneracies are on logarithmic scale for 
clarity. 

Fig. 3 Standard deviation of the mean field as a function of total number of sites for 
the CML with function f(x) = fix(l — x), with /i = 4.0 and e = 0.1. Similar 
behaviour also obtains for other values of /i, e. 

Fig. 4 Probability distributions of the CML at a central point and at the boundary 
for the same set of parameters as in Fig. 3. The lattice has k = 6, and the first 
10 4 iterates are treated as transients. The distribution is clearly very different 
for these two sites. 



21 



Fig. 1 




22 



Fig. 2 
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Fig.3 
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Fig.4 
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